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A multidimensional cosmological model describing the dynamics of n + 1 Ricci-flat factor-spaces Mi in the presence 
of a one-component anisotropic fluid is considered. The pressures in all spaces are proportional to the density: 
Pi = Wip, i — 0, ...,n. Solutions with accelerated expansion of our 3-space Mo and small enough variation of the 
gravitational constant G are found. These solutions exist for two branches of the parameter wo ■ The first branch 
describes superstiff matter with wo > 1, the second one may contain phantom matter with wo < —1, e.g., when G 
grows with time. 



1. Introduction 

There are many hot spots in the gravitational interac- 
tion 1 . Among them one may point out such problems 
as acceleration of the Universe expansion, description 
and detection of strong field objects (black holes, worm- 
holes etc.) and gravitational waves, near-zone experi- 
ments, such as equivalence principle tests, second-order 
tests, rotation and torsion effects of general relativity 
etc. Within the last block, we want to stress a special 
role of experiments to measure the value of the grav- 
itational (or Einstein) constant and its possible varia- 
tions. These experiments already belong to new gener- 
ation ones since they are testing not only the gravita- 
tional interaction but also some predictions of unified 
models and theories. A special role in these activities is 
played by space experiments, and this role will increase 
in the future. Modern cosmology and its observational 
part already became an arena for testing predictions of 
high energy physics. All this leads to a fundamental 
role of gravity in the present investigations, it is still 
a missing link in unified theories. Fundamental physi- 
cal constants, relations between them and their possible 
variations are a reflection of the situation with unifica- 
tion [lX2uA\- 

We will be mainly interested in the gravitational con- 
stant and its possible variations. There are three prob- 
lems connected with the Newtonian gravitational con- 
stant G d]: 

1. Absolute value of G. 

2. Possible time variations of G . 

3. Possible range variations of G, or new interac- 
tions (forces). 

The oldest one is the problem of possible tempo- 
ral variation of G , which arose due to papers of Milne 
(1935) and Dirac (1937). In Russia, these ideas were de- 
veloped in the 60s and 70s by K.P. Staniukovich [51l2j, 
who was the first to consider simultaneous variations of 
several fundamental constants. 
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Our first calculations based on general relativity with 
a perfect fluid and a conformal scalar fleld [5] gave G/G 
at the level of 10^-'^^ — lO^^'^ per year. Our calculations 
in string-like [7j and multidimensional models with per- 
fect fluid f8] gave the level 10~^^ , those based on a gen- 
eral class of scalar-tensor theories [11] and simple mul- 
tidimensional model with p-branes [inilll] gave for the 
present values of cosmological parameters 10"^'^ — lO"^** 
and 10~^^ per year, respectively. Similar estimations 
were made by Miyazaki within Machian theories [13j 
giving for G/G the estimate 10~^^ per year and by Fu- 
jii — on the level 10~^^— 10~^^ per year [Uj. Analysis of 
one more multidimensional model with two curvatures 
in different factor spaces gave an estimate on the level 
10^^^ [15 . Here we continue our studies on variation of 
G in another multidimensional cosmological model. 

2. The model 

We consider a cosmological model describing the dy- 
namics of n Ricci-fiat spaces in the presence of a 1- 
component "perfect-fiuid" matter [l6j. The metric of 
the model 



g = - exp[2j{t)]dt (g) dt ^ 

is defined on the manifold 
Af = R X A/o X . . . X Mn 



E 

4=0 



e^p[2x\t)]g' 



(2.1) 



(2.2) 



where Mi with the metric is a Ricci-fiat space of di- 
mension di, 1 = 0,. ..,n; n > 2. The multidimensional 
Hilbcrt-Einstein equations have the form 



Rn 



(2.3) 



where is the gravitational constant, and the energy- 
momentum tensor is adopted as 

{T'J) = diag(-p, poC' • ■ • ' P-'^D: (2-4) 
describing, in general, an anisotropic fluid. 
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We put pressures of this "perfect" fluid in all spaces 
to be proportional to the density, 



Piit) = (1 - Ui/d^)p{t), 



(2.5) 



where Ui = const , i = 0, . . . ,n. We also put p > . 

We impose also the following restriction on the vec- 
tor u = (ui) e : 



(m, u)* < 0. 

Here, the bilinear form 
defined by the relation 

(u, u)* — G^-'uiVj, 

u,v £ M"+^ , where 

S'^ 1 



2-D 



(2.6) 
M is 

(2.7) 
(2.8) 



are components of the matrix inverse to the matrix of 
the minisuperspace metric pO[[2T| 

Gij — diSij — didj. (2.9) 

In (j2.8p . D = 1 + X]"=o '^i ^^'^ total dimension of the 
manifold M (g^)- 

The restriction (12.61) reads 



1 / " \^ 

-J,[J2^^) <0. (2.10) 



2-D 



3. Solutions with power-law scale 
factors 

Here, we consider a special family of "power-law" solu- 
tions from |161ll7j with the metric written in the syn- 
chronous time parametrization 
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-dts (g) dts + ^ai{t,)g\ 



(3.1) 



i=0 



Solutions with a power-law behaviour of the scale factors 
take place for 



Here and below the vector 



(3.2) 



(3.3) 



corresponds to the A-term fluid with pi — —p (vacu- 
umlike matter). 

In this case, the solutions are determined by the met- 
ric (13.11) with the scale factors 



(3.4) 
(3.5) 
(3.6) 



and the density 

— 2( u, u )* 



-u,M)2^2• 



Here 



where = G^^Ui and Ai are positive constants, i = 
0,...,n. 

The model under consideration was integrated in [16j 
for (u,it)* < 0. The solutions from 16J were gener- 
alized in [17j to the case when a massless minimally 
coupled scalar held was added. Families of exceptional 
solutions with power-law and exponential behaviours of 
the scale factors in terms of synchronous time were sin- 
gled out in [17 1 and correspond to a constant value of 
the scalar field: (p = const. When the scalar field is 
omitted, we are led to solutions presented above (in T6j 
these solutions were originally written in the harmonic 
time parametrization). It may be verified that the ex- 
ceptional solutions with power-law dependence of scale 
factors are also valid when the restriction p.6p is omit- 
ted. Moreover, it may be shown that for ( u, m ), =0 the 
power-law solutions coincide with the vacuum Kasner- 
like solution from [22]. In this case, the matter source 
vanishes since p = in (|3.5[) . 

4. Acceleration and variation of G 

In this section, the metric is assumed to be fiat, 
and do = 3. The subspace (Mo,g°) describes "our" 
3-dimensional space and {Mi,g^) internal factor-spaces. 

We are interested in solutions with accelerated ex- 
pansion of our space and small enough variations of the 
gravitational constant obeying the present experimental 
constraints, see 1101: 



\G/{GH)\it,o) < 0.1, 
where 



ao 



(4.1) 



(4.2) 



is the Hubble parameter. We suppose that the internal 
spaces are compact. Hence our 4-dimensional constant 
is (see |H|) 



G = const •]^"_^(ar''' 



(4.3) 



We will use the following explicit formulae for the 
contravariant components: 

«^ = G^^^. = ^ + ^E".- (4.4) 

j=o 

and the scalar product reads 

^ d^ D - 2^ D ~2^^ ' ^ ' 

1=0 j=0 1=0 



4.1. Power- law expansion with acceleration 

For solutions with power-law expansion, an accelerated 
expansion of our space takes place for 



> 1. 



(4.6) 
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For D ~ A, when internal spaces are absent, we get 
ly" = 2/{6-uo), (4.7) 

(4.8) 
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(u^^^ - M, u), = |(wo - 6)uo 7^ 0, 

which impHes uq ^ and uq ^ 6 (here (w, u)* = 
— < 0). The condition j/'^ > 1 is equivalent to 
4 < uo < 6 , or, equivalently, 

-p<p<-p/3, (4.9) 

which agrees with the well-known result for = 4. (We 
note that recently special 5-dimensional power-law solu- 
tions (e.g., with acceleration) were considered in |23j). 
For power-law solutions we get 



G EU ^'d. 



G ts 
and hence 

G/{GH) = -^ 

1 A' 



do ts 



y'd, EE S. 



(4.10) 



(4.11) 



The constant parameter S describes variation of the 
gravitational constant and, according to (|4.ip , 



|(5|<0.1. (4.12) 
It follows from the definition of i^* in (13.611 that 



1 



or, in terms of covariant components (see (j4.4p ) 
(i?-4)uo-2Er=i«. 



5 = - 



i(5 - d)uq + x;r=i 



(4.13) 



(4.14) 



Thus the relations (gjH), (|iT^ . (|iTi)) and the 

constraint p.2p determine a set of parameters Ui com- 
patible with the acceleration and tests on G-dot. 

In what follows we will show that these relations do 
really determine a non-empty set of parameters Ui de- 
scribing the equations of state. 

4.1.1. The case of constant G 

Consider the most important case (5 = 0, i.e., when the 
variation of G is absent: G = Q. 

Indeed, there is a tendency of lowering the upper 
bound on G . Moreover, according to arguments of T9], 
5 < 10^"* . This severe constraint just follows from the 
identity 



G/G = a/a 



(4.15) 



that takes place in multidimensional models. Here a is 
the fine structure constant. 

Isotropic case. First we consider the isotropic case 
when the pressures coincide in all internal spaces. This 
takes place when 



1, 



(4.16) 



For pressures in internal spaces we get from (j2.5p 

p, = {l-v)p, z = l,...,n. (4.17) 

Then we get from ((2T0| and (|4?5)) 
1 



\{d - l)uQ + 2duov - 2dv% (4.18) 



1 



\2uo + 2dv 



2~D' 

+ \{d-l)ul-2duov + 2dv^]. (4.19) 



Here and in what follows we denote d = D ~ 4. 
For 5 = , we get in the isotropic case 



(4.20) 



V = uo/2, 

or, in terms of pressures, 

= (3po - p)/2, i = l,...,n. (4.21) 

Substituting ((420)) into ((4T8)) and ((4T9)) we get 

{u,u), = ~ul/6, (4.22) 
(u(^) ='Uo(uo-6)/6. (4.23) 

Remarkably, we obtain the same relations as in Z? = 
4 case (see the remark above). For our solution, we 
should put Up ^ and Uq ^ 6. 

Using (jil^ and (|i?jDl) we get m° = ~uo/6 and 

= for z > 0, hence i^i = for i = 1, . . . , rt, i.e., all 
internal spaces are static. 

The metric p.ip reads in our case 



-dts ® dts + Altl""g^ + ^ Ajg\ 



(4.24) 



where Ai are positive constants, and 

I'" = 2/(6- uo). (4.25) 

We see that the power is the same as in case D = 4. 
For the density we get from (|3.5p 



p 



12 



("0 - 6)2^2 ■ 



(4.26) 



Thus the equations of state (|2.5[) . with relations 
(jiTTO)) and imposed, lead to the solution 

(|4.26p with a Ricci-fiat (e.g., flat) 3-metric and n static 
internal Ricci-fiat spaces. For 



4 < uq < 6 



(4.27) 



or, equivalently, —p < pq < — p/3, we get an accelerated 
expansion of "our" 3-dimensional Ricci-fiat space. 

Aisotropic case. Consider the anisotropic (w.r.t. in- 
ternal spaces) case with i5 = 0, or, equivalently (see 

631) , 



{D-4)uo = 2j2'' 



(4.28) 
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This implies 



where 



u, u), = huo{uo - 6) 



1^0 



A, 



n 2 1 / " 



>o, 



(4.29) 
(4.30) 



(4.31) 



The inequahty in (|4.31[) can be readily proved using 
the well-known Cauchy-Schwarz inequality: 



> 



(4.32) 



Indeed, substituting bi = ^/di and Ci = Ui/^/di into 
(|4?32)) . we get (|4?3T1) . The equality in ((432)) takes place 
only when the vectors (bi) and (q) are linearly depen- 
dent, that for our choice reads: Ui/^fdi — v\fdi where 
V is constant. Thus A = only in the isotropic case 
(|4.16p . In the anisotropic case we get A > . 
In what follows we will use the relation 



where 



= 3 ± \/9" 



6A 



(4.33) 



(4.34) 



are roots of the quadratic trinomial (|4.29p obeying 
W(7 < 0, u+ >6 for A > 0. (4.35) 



It follows from ((i?^ that 
2uo 



-uq/6 and hence 



uo{uq — 6) — 6A 



(4.36) 



(here uq ^ u^). 

The function iy'^{uo) is monotonically increasing: (i) 
from to -|-oo in the range (— oo, Uq) ; (ii) from — oo to 
-|-oo in the range {uq,Uq); (iii) from — oo to in the 
range {uq,+oo). This behaviour of j/°(uo) in each of 
the three ranges simply follows from the relation 



2{ul + 6A) 



duQ (uo - u^)'^ (uo - Uq )^ 



> 0. 



(4.37) 



The accelerated expansion of our space takes place 
when I/" > 1 , or, equivalently, when either 



(A) Uo e (mo*>Wo )' 

(B) uoeiu^^,u^), 

where 



2± V4T6A. 



In terms of the parameter wq , 
Po = Wo/5, Wo = 1 - wo/3. 



(4.38) 
(4.39) 



(4.40) 
(4.41) 



these two branches read: 
(A) 



Wo =^l + |A<wo< i + l^l-f |A = Wo,, (4.42) 
(B) 



Wo+ = -^1 + |A < wo < i - 1^1 + |A = w+. (4.43) 

The first branch (A) describes superstiff matter (wo > 
1 ) with negative density. Indeed, p < follows from 
(P31) and > 0, see (4.30). 

The second branch (B) corresponds to matter with 
a broken weak energy condition (since wo < — | ) and 
positive density (since < 0). This matter is 

phantom (i.e., wo < — 1) when A > 2. For A < 2 the 
interval {wq,Wq^) contains both "phantom" (wo < — 1) 
and "non-phantom" points ( wo > ~ 1 ) • 

4.1.2. The case of varying G 

Now we consider another important case S ^ 0, i.e., 
when G 0. In what follows, we use the observational 
bound (|4.12p : \5\ < 0.1, stating the smallness of S. 
Using (|4.14p . we get 

n ^ 

= -d6Mo, (4.44) 

1=1 

where d = D — 4 and 

For the scalar product we get from (|4.44p 



{u^^^ -u,u)^::^ lAul~ Buo- A, (4.46) 
{u,u), = -lAul + A, (4.47) 

where A was defined in (|43T|) (see (f2T0|) and ([431)). 



A 



1 d,V 1 



B = 



1 



d + 2 



(2 + db). 



Using (|4.45p . we obtain the explicit formulae 
(d + 2)(52 



A = A{5) = 1 
B = B{5) 



12d(l- (5/2)2^ 
1 -(5/3 



1-5/2' 



(4.48) 
(4.49) 

(4.50) 
(4.51) 



It should be noted that, due to \5\ < 0.1, A is positive, 
A > 0, and close to unity: |A- 1| < |10"2. 

For the contravariant component vP we get from 
(113) and (H^ii)) : 



u" = -Cuo/6, 
where 

G = G{S) =3B-2 = 1/(1 - 5/2). 



(4.52) 
(4.53) 
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It follows from (11:17)) and that (see HjJ^ ) 

2Cmo 

~ Aul/6~ Buq- A' 

Here uq ^ where 



(4.54) 



4 = = ^(3B± v/9B2T6AA) (4.55) 



are roots of quadratic trinomial (|4.46p . 
In what follows we will use the identity 



Aul - 4uo - 6A 
Au\ - 6Buo - 6A ■ 



(4.56) 



Isotropic case. Let us consider the isotropic case 
(|4?T6)) . Then we obtain from (1444)) 



V = dbuo/2- 
or, in terms of pressures 
1, 



(4.57) 



p,; = -[36po + (2-36)p], i = l,...,n. (4.58) 



For scalar products we get 

{u^^^ - u,u)^ = uo{Auq - 6B)/6. 



(4.59) 
(4.60) 



For our solution, wc should put uq ^ and uq ^ 
QB/A. The metric p.ip reads in our case 

n 

g = -dts ® dts + Altf'g° + tf (4-61) 



where Ai are positive constants. 



2C 



Alto - 6S ' 



and 
25 



d{\ - (5/2)(Awo - 65)' 
Tl 
uq5 



(4.62) 
(4.63) 



i = 1 , . . . , n . The last formula follows from (|3.6p and 



6d(l - 5/2)' 



(4.64) 



We see that the power does not coincide, for (5 7^ 0, 
with that in case D — A. 

For the density, since A > , we get from ()3.5p 



12A 



> 0. 



(4.65) 



The accelerated expansion condition for our 3D 
space, i^" > 1 , reads in this case 



4 6B(S) 
<ua< 



AiS) ^ - ^ AiS) 
or, equivalently, in terms of wq ()4.41|) (po = "u^op) 



(4.66) 



3A(5) 



= i/;+((5).(4.67) 



For 5 > 0, we get an isotropic contraction of the 
whole internal space Mi x . . . x M„ . In this case, 
Wq (6) < — 1 , and hence phantom matter may occur 
with the equation of state close to the vacuum one since 

5(1 + S/d) 



w+i5) + l = - 



3[l-5 + {d-l)d^/{6d)y 



(4.68) 



For small 6 we have w^{S) = — 1 — 6/3 + 0{6'^) . 

For S < we get an isotropic expansion of the whole 
internal space. Then 'Wq{S) > — 1, and phantom mat- 
ter does not occur. In both cases Wq^{S) < — 1/3 and 
wl{S) + 1/3 = 0{S^). 

Anisotropic case. Now we consider the anisotropic 
case A > when 5 7^ 0. Using ()4.56p . we obtain 



{uq - U^)(uo - Uq) 

where ~ Uq{6) were defined in (4.55) and 



(4.69) 



"o"; = 4M = 2 ± ^4 + 6^(5) A. (4.70) 

Accelerated expansion of our 3-dimensional space 
takes place when i/'^ > 1 , or, equivalently, when either 

(A) uoe{uo,i5),UoiS)), or (4.71) 

(B) uoe{ul{S),u+iS)). (4.72) 

In terms of the parameter wq (pa = Wop, Wq = 1 — ^) 
these two branches read: 

(A) w^{S) <wo <w^^{S), 

(B) w+{S)<wo<wliS), 
where 

w^iS) = 1^4{S)/3, 
wiid) = l-4*W/3. 
For small 6 we have 

6V V9 



6A 



(4.73) 
(4.74) 

(4.75) 
(4.76) 

0(52), (4.77) 



(4.78) 

Thus for small S the lower and upper bounds on wq 
have a small deviation from those obtained for 5 = 0. 
For small 5, the upper bounds shift only by 0{S^) , while 
the lower bounds shift by 0(5) . 

The first branch (A) describes superstiff matter 
Wq > 1 since Wq (6) > 1 due to Ug{5) < 0. It may 
be shown that the density is negative in this case since 
> 0. 

For branch (B) we get for the upper bound 
«;+(5)<-l/3 (4.79) 

due to Uq^{S) > 4. For the lower bound we find that 
w+(5) < -1 (4.80) 

only if 

A > 6[A(5) - B{S)] = -5/(1 - 5/2)2. 

This is a condition on the appearance of "phantom" 
matter. For 5 > this inequality is valid, but for 5 < 
it is satisfied only for a big enough value of anisotropy 
parameter A, see ()4.8ip . 
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5. Conclusions 

We have considered multidimensional cosmological mod- 
els describing the dynamics of n + 1 Ricci-flat fac- 
tor spaces Mi in the presence of a one-component 
anisotropic fluid with pressures in all spaces propor- 
tional to the density: pi ~ Wip, i = 0, ...,n. Solutions 
with accelerated expansion of our 3-dimensional space 
Mq and small enough variation of the gravitational 
constant G were found. These solutions exist for two 
branches of the parameter wq : (A) from (|4.73p and (B) 
from (4.74). Branch (A) describes superstiff matter with 
wq > 1 while branch (B) may contain phantom matter 
with wq < — 1 . The second branch obviously contains 
phantom matter states when (i) 6 > 0, i.e., G increases 
with time, or (ii) 6 = and A > 0, i.e., if G is constant 
and the expansion (or contraction) of the factor spaces 
is anisotropic. If (iii) (5 < 0, i.e., G decreases with time, 
phantom matter appears if the anisotropy parameter A 
is big enough (see (|4.8ip .) 
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